Time-optimal CNOT between indirectly coupled qubits in a linear Ising chain 
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We give analytical solutions for the time-optimal synthesis of entangling gates between indirectly 
coupled qubits 1 and 3 in a linear spin chain of three qubits subject to an Ising Hamiltonian 
interaction with a symmetric coupling J plus a local magnetic field acting on the intermediate 
qubit. The energy available is fixed, but we relax the standard assumption of instantaneous unitary 
operations acting on single qubits. The time required for performing an entangling gate which 
is equivalent, modulo local unitary operations, to the CNOT(l,3) between the indirectly coupled 
qubits 1 and 3 is T = ^/3/2 J -1 , i.e. faster than a previous estimate based on a similar Hamiltonian 
and the assumption of local unitaries with zero time cost. Furthermore, performing a simple Walsh- 
Hadamard rotation in the Hlibert space of qubit 3 shows that the time-optimal synthesis of the 
CNOT ± (l,3) (which acts as the identity when the control qubit 1 is in the state |0), while if the 
control qubit is in the state |1) the target qubit 3 is flipped as |±) — > |=p)) also requires the same 
time T. 
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I. INTRODUCTION 

Quantum optimal control theory is by now a very well 
studied subject, both theoretically and experimentally, 
with increasing applications in the field of quantum com- 
puting and information (for an updated review see, e.g., 
|l| and references therein). Quantum optimal control 
techniques aim, e.g., at finding either the best quantum 
state evolution or the best unitary operator evolution 
with respect to some fixed cost, which can be assumed 
to be a fidelity with respect to a target (a quantum state 
or a gate), the purity of the target state, etc. Time- 
optimal quantum computation, where the cost to be op- 
timized is the time to achieve a given quantum evolution, 
has also recently become a hot topic in quantum con- 
trol and quantum information theory j^-[23|. For exam- 
ple, minimization of the physical time to achieve a given 
unitary transformation is relevant for the design of fast 
elementary gates and provides a more physical ground 
to describe the complexity of quantum algorithms than 
the standard concept of gate complexity, which gives the 
number of elementary gates used in a quantum circuit 
[23j . Also the literature in time-optimal quantum control 
is rapidly growing and it is not the purpose of this paper 
to give a complete review. For instance, Refs. 0-0 dis- 



cuss the time optimal generation of unitary operations 
for a small number of qubits using Lie group methods, 
the theory of sub Riemannian geometry, the Pontryagin 
maximum principle and assuming that one-qubit opera- 
tions can be performed arbitrarily fast. The time-optimal 
synthesis of unitary transformations (quantum gates) be- 
tween two coupled qubits has been thoroughly discussed 
0-0 ■ The time-optimal evolution of quantum states 
for qubits and qutrits has been also investigated, e.g., in 
[Iol|-jllj|. that of a 2-level dissipative system, e.g., in fl2| . 
while the time-optimal generation of cluster states has 
been considered by [l3j]. The extension of the Lie alge- 
braic methods to the coupling of slow and fast systems 
can be found in [bfj ]. Earlier bounds on the time com- 
plexity of generating 2-qubit unitary gates can be found 
in |l5| , while lower bounds on the time complexity of n- 
qubit gates are given in [l6| and upper bounds on the 
time complexity of certain n-qubit gates via several cou- 
pling topologies are numerically described in [TtJ ■ The re- 
lationship between time complexity and gate complexity 
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has been also investigated in [18j. Nielsen et al. [19| 
proposed a criterion for optimal quantum computation in 
terms of a certain geometry in Hamiltonian space, and 
showed that the quantum gate complexity is related to 
optimal control cost problems. An adiabatic solution to 
the optimal control problem in holonomic quantum com- 
putation has also been found in (22[. Finally, numerical 
methods for the design of optimal quantum control evo- 
lutions have been proposed, based on the gradient ascent 
algorithm p4| and on the non linear Krotov algorithm 

Most of the above works share the assumption that 
one-qubit gates have zero time cost. The present au- 
thors described a theoretical framework for time-optimal 
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quantum computing based on the action principle where 
such an assumption is not necessary, and named it the 
quantum brachistochronc (from hereon abbreviated as 
QB (26[). The variational principle is formulated for 
the time-optimal evolution of a quantum system whose 
Hamiltonian is subject to a set of constraints (e.g., a finite 
energy or magnetic field, certain qubit interactions are 
forbidden) and defines a boundary value problem with 
fixed initial and final quantum states (or unitary trans- 
formations) . The QB has been studied for quantum state 
evolution in the case of pure [27j and mixed states f28j . 
and for the optimal realization of unitary transformations 
between the identity and a given target quantum gate 
[H| . The latter is particularly relevant for the standard 
quantum computation paradigm since a whole algorithm 
may be reduced to a sequence of unitary transformations 
between intermediate states and a final measurement. 
The more realistic situation where the target quantum 
state (gate) can be reached within a finite, tolerable er- 
ror (a fidelity larg er than a specified value) has also been 
addressed in [301 ] . The quantum brachistochrone prob- 
lem always reduces to solving a fundamental equation, 
which can be easily written down once the constraints for 
the Hamiltonian of the quantum system are known, with 
given initial boundary conditions, and an equation for the 
Lagrange multiplier which enforces the dynamical law for 
the quantum system (the Schrodinger equation for closed 
systems or, e.g., a master equation for Markovian open 
systems), with given final boundary conditions. Two of 
us also studied numerically the time complexity of gen- 
erating unitarics acting on n qubits via a Hamiltonian 
which contains only one and two qubit interaction terms 
[3l| . Our research also triggered an amount of related 
works [HI- [i^. For example, the authors of [33j]- 3f| con- 
sidered the problem of the generation of multipartite en- 
tanglement during the QB evolution of quantum states 
and unitaries, while those of (36|- [4(| studied the QB in 
the context of non Hcrmitian quantum mechanics (for a 
review of the latter works see, e.g., (4l|V 

More recently, Khaneja et al. [42[ considered the prob- 
lem of the efficient synthesis of the controllcd-NOT gate 
(CNOT(l,3)) between two qubits indirectly coupled via 
an Ising type coupling to a third qubit, and where the 
single qubits can be separately addressed via instanta- 
neous local unitaries. Their implementation requires a 
time T ~ 1.253J -1 ~ \p^T^ J 1 -, where J is the Ising 
coupling between the qubits (1, 2) and (2, 3) in a linear 
coupling topology, and it was related to the computation 
of a geodesic on the surface of a sphere with a special met- 
ric and to the synthesis of a particular entangling gate 
called U( 3 . Subsequently, similar methods were used to 
extend these results to the case of unequal Ising couplings 
between the indirectly coupled qubits 43| and to the case 
of a spin chain with n qubits [4J] . The time-optimal syn- 
thesis of interactions between qubits which are indirectly 
coupled via an intermediate qubit is a typical scenario in 
quantum information processing (e.g., distant spins on 
spin chains (45j . bus-qubit electrons mediating interac- 



tions of two nuclear spins [46[ or NMR experiments cor- 
relating the frequencies of indirectly coupled qubits [I?} ) . 
The synthesis of such indirect couplings is usually time 
costly (e.g., it may require concatenation of two-qubit 
operations on the directly coupled qubits) and therefore 
it is also easily prone to decoherence effects and a degra- 
dation of the gate fidelity, both critical effects in practical 
implementations of quantum computing. 

In this work we further investigate the problem of the 
time-optimal generation of the same entangling quantum 
gates, U( 3 and CNOT(l,3), between two indirectly cou- 
pled qubits in three-qubit systems with linear coupling 
topologyand with similar Ising coupling Hamiltonian as 
that of [42|. We find that, if the available Hamiltonian is 
made of an interacting piece of the Ising type exactly as in 
the model of [42[ plus a local magnetic field acting on the 
intermediate qubit, and is subject to the constraint of a 
finite energy, then the C/f 3 gate can be optimally realized 
in a time T = J -1 , which is faster than the time 

required via the construction of [42j . Furthermore, using 
a slightly modified interaction Hamiltonian (obtained via 
a simple change of basis for the Hilbcrt space of one of the 
indirectly coupled qubits, qubit 3) and the same energy 
constraint, we find that the gate CNOT ± (l,3) (which 
acts as the identity when the control qubit 1 is in the 
state |0), while if the control qubit is in the state |1) the 
target qubit 3 is flipped as |±) — > |=p)) can be also opti- 
mally generated in the same time T = \/3/2 J -1 as the 
entangling gate t/f 3 . The paper is organized as follows. 
In Section I we briefly review the main features of the 
QB formalism for the time-optimal synthesis of unitary 
quantum gates. In Section II we discuss the problem of 
the efficient generation of the gate U{ 3 between the in- 
directly coupled boundary qubits of a three-linear qubit 
system subject to an Ising interaction and a local con- 
trol available for the intermediate qubit, when a finite 
energy is available. Section III is devoted to the study 
of the time-efficient generation of the CNOT ± (l,3) gate 
with the slightly modified interaction Hamiltonian plus 
the same local operation on qubit 2 and the same energy 
available. Finally, Section IV is devoted to the summary 
and discussion of our results. 



II. QUANTUM BRACHISTOCHRONE 

We want to find the time-optimal way to generate a 
target unitary operation Uf £ U (N) (modulo physically 
irrelevant overall phases) by controlling an Hamiltonian 
H{t) and evolving a unitary operator U(t), both obeying 
the Schrodinger equation. We assume that H is control- 
lable within a certain available set, dictated either by the- 
oretical conditions (e.g., only certain interactions among 
qubits are allowed) or by experimental requirements (e.g., 
a finite energy or a finite magnetic field). At least the 
'magnitude' of the Hamiltonian must be bounded, other- 
wise any gate U / might be realized in an arbitrarily short 
time simply by rcscaling the Hamiltonian [29[ . Physically 
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this corresponds to the fact that one can afford only a fi- 
nite energy in the experiment. 

The tim e-op timality problem is formulated using the 
action [H, : 

S{U,H;a,A,Xj) := [ dT[a + L s +L c ] (1) 
Jo 

(A,i%Ui-aH), (2) 

aJ2^f j (H), (3) 



Lc 



where (A, B) := Tr(A'B) and the Hermitian operator 
A(t) and the real functions Xj(r) are Lagrange multipli- 
ers |48j |. The quantity a is the time cost, and it may 
be interpreted [28j as a positive independent dynamical 
variable (a " lapse" function) which measures the physical 
time t := J a(r)dr lapsed in each infinitesimal interval 
dr of the parameter time r. Variation of Ls by A gives 
the Schrodinger equation: 



or 



U(t)=Te- i fi adt , 



where T is the time ordered product. Variation of Lc by 
A j leads to the constraints for H: 



The quantum brachistochrone together with the con- 
straints define a boundary-value problem for the evo- 
lution of the unitary operator U (t) with fixed initial 
(U(t = 0) = f, where 1 is the identity matrix) and fi- 
nal conditions (U(t = T) = Uf, where T is the optimal 
time duration necessary to achieve the target gate Uf). 
The quantum brachistochrone is a set of first-order (non 
linear) differential equations which can always be solved 
in principle, e.g. numerically, and it is universal, as it 
holds also in the case of time-optimal evolution of pure 
[27[ and mixed [28| quantum states. 

More in details, for a given target gate Uf, the pro- 
cedure to find the optimal Hamiltonian H and the op- 
timal time duration T consists of the following stages: 
(i) specify the constraint functions fj(H) for the avail- 
able Hamiltonian; (ii) write down and solve the quantum 
brachistochrone ([9]) together with the constraints ((SJ) to 
obtain Hoprit); (hi) integrate the Schrodinger equation 
05) with U(0) = 1 to get U pr{t); (iv) fix the integra- 
tion constants in Ho pt (t) by imposing the condition that 
Uopt(T) equals Uf modulo a global (physically irrele- 
vant) phase, i.e., 



Uopt(T) = 
where \ is some real number. 



Uf, 



(10) 



f J (H)=Q. 



(5) 



In particular, the finite energy condition for a system of 
log V qubits can be written as: 



fo(H) := ±[Tr(ff 2 



Nco A 



0, 



(6) 



where u is a constant. Since overall phases are irrele- 
vant in quantum mechanics, it is natural to consider the 
time-optimal evolution of unitary operators belonging to 
the group U{N)J U(l) ~ SU(N) (i.e., the dynamics is 
generated by a traceless Hermitian Hamiltonian) . 
Then, we introduce the operator: 



F 



dLr 



dH ' 

and from © and (J7J we obtain 

F = \ H + F'. 



(7) 



(8) 



If the other constraints are linear and homogeneous in 
H, i.e. F' = J2j Ajffj with gj e su(N), then we have 
Tt(HF') = and it is easy to show [H| that the 
Lagrange multiplier Ao in (|8]) is a constant. 

From the variation of S by a, upon using (JB]) and ([5]). 
one gets 1 = Tr(HF) = \ Tr(H 2 ) = A iVw 2 , which 
determines the constant Ao- 

Finally, variation of S by U and some elementary al- 
gebra give the quantum brachistochrone equation 



III. ISING HAMILTONIAN AND 
TIME-OPTIMAL ENTANGLER Uf 3 

Now we apply the general QB formalism summarized 
in the previous Section to the case of a physical sys- 
tem of three qubits represented by three spins (labeled 
by a superscript a £ {1,2,3}) interacting via an Ising 
Hamiltonian with time independent couplings J12 , J23 
and subject to a local and controllable magnetic field 
Bf(t) (i = x, y, z). In other words, we choose the three- 
qubit Ising Hamiltonian, 



MtKo-t] +B(t)- a 2 , (11) 



where we have used the simplified notation, e.g. 
o~i ® o-j ® 1, ofer;- := 1 ® < jj ® o-j, 0} :— 1 ® crj ® 1 and 
crj are the Pauli operators [2j]. We further assume that 
the Ising couplings in (fTTj) are equal and (a positive) con- 
stant, i.e. J 12 {t) = J 23 (t) := J > (HI. This is formally 
enforced via the Lagrange multipliers v zz (t) and p zz (t) 
and, respectively, the associated constraints 



h := MHa 2 z a 3 z ) 



AttJ 
4ttJ 



Moreover, the finite energy condition ^ reads 



B =ur — 



= const. 



(12) 
(13) 



(14) 



(9) 



We note that the interaction part of the Hamiltonian pT|) 
is exactly the same as H c in eq. (2) of (42[, and that the 
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local term in (fTTj) corresponds to the terms Ha and Hb 
of [42[ . Finally, the form (fTTj) of the physical Hamiltonian 
is guaranteed by the operator 



F'(t) = y^\, M l^;^l 

i,j,k 

+ Y^Mt)oiOi + + fHj(t)<>iOj] 

+ E[* 1+ 6(')°f] ( 15 ) 

where Ayfc(i), My(t), (i) , Tfc (i) and are La- 

grange multipliers, and the indices € {x,y,z}. 

Our task is then to solve the quantum brachistochronc 
equation ©. Comparing the coefficients of the genera- 
tors of su(8) on both sides of ©, with F given by (J5J, 
7? given by (jlip and F' given by (fT5j) , we find that the 
relevant quantum brachistochronc equations are: 

B x = -TrJ(v zy + p yz ) 
B y = ttJ(v zx + p xz ) 
B z = 



Vz 



= -InJX 



zyz 



2{B Z V Z y - ByV ZZ )\ 



v zy = ~kJ\ zxz + 2(B z v zx — B x v zz ) 

v zz = 2(B x v zy — B y v zx ) 

Pxz = -[nJ^zyz + 2(B z p yz - Byp zz j\ 

Pyz = kJ\ zxz + 2(B z p xz - B x p zz ) 

pzz = 2(B x p yz — B y p xz ) 

K.xz = —irJ(v Z y + p yz ) + 2(B y \ zzz - B z \ zyz ) 



kJ (v zx + p xz ) - 2(B X \ Z 



X ZZZ — 2{B X \ Z y Z — ByX ZXZ ) 



(16) 



As recalled in the previous Section, the first step in ad- 
dressing the QB problem is to solve the fundamental 
equation ([9]) for the time-optimal Hamiltonian Hopr(t)- 
In other words, we need to solve cqs. (|16[) together with 
the constraints (Ti"2"T) - (|14[) and with fixed parameters J,lu. 
From the integral of the motion B z = const, eqs. ([To]) 
and the energy constraint (|14[) we immediately obtain 
that B% + B^ := Bq = const while, from eqs. ([To]) , we 
further find the integrals of the motion v zz + p zz = const, 
\ zzz = const. Exploiting these integrals of the motion 
and after some lengthy but elementary algebra, we find 
that the general and non trivial solution of (|16|) is given 
by the Hamiltonian (fTTj with the time-optimal magnetic 
field 



B Q cos (t) \ 
B Q sm0(t) j 



(17) 



precessing around the z-axis with the frequency f2, where 
0(t) := Qt + 0(0) and f2 and 0(0) are integration con- 
stants. 

The next step is then to integrate the Schrodinger eq. 
([JJ for Uopr(t), given that Hopr(t) is expressed by eqs. 



(ITTl) and ([IT]). For this purpose, we exploit the following 
well known property for the rotation of the Pauli matri- 
ces: 



e % 2 a 'a x e l 2 a * = cos 6(t) a x + sin 0(t)a y , 
and we rewrite the time-optimal Hamiltonian as 

He 

where we have introduced the constant operator 
"ttJ. 



■Sit) 2 -Bit) 2 



Hn := Bno„, 



-(al+al) + B z 



(18) 



(19) 



(20) 



Furthermore, defining the transformed unitary operator 
U(t):=e i!S P a iu(t), (21) 



we easily check that, since U(t) should obey the 
Schrodinger equation (j4]), U(t) should also satisfy 



if=HU 
dt 

with the time-independent Hamiltonian 

n n 

H := Ho - —o z = const. 



(22) 



(23) 



We note that the constant Hamiltonian H is diagonal in 
the 1,3 qubit subspace, i.e. 



H = B ai + B^a 2 z = const, 



(24) 



where we have introduced the operator (in the 1,3 qubit 
subspace) 



R 13 



-(n-2B z )l + nJ Diag[l, 0,0,-1]. (25) 



Then, solving eq. together with cq. (gU) for U(t) 
and finally inverting (|2"Tj) . it is easy to check that the 
time-optimal unitary operator UoPT(t) evolves as 



J7opt(*) 



e 2 



(26) 



In particular, the exponential of the constant Hamilto- 
nian appearing on the right hand side of eq. (|2l)|) is also 
diagonalized (in the 1, 3 qubit subspace) and can be ex- 
panded as 



= C 1 D 3 (t)~iS 1 D 3 (t)H, 



(27) 



where we have introduced the following operators (acting 
in the 1,3 qubit subspace) 

S£(t) := Diag[8 + (t),so(t),so(t),s-(t),] (28) 
C^(t) := Diag[c+(t),co(*),Co(i),c_(*)], (29) 
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which depend upon the functions 
sin u±t 



s±(t) := 
c±(t) := 
and the constants 



s (t) :-- 



sin wq t 



cosw±i; co(t) := coswoi, 



:= Bf, 



Bl 



[SI - 2(B Z ± ir J)f 
(Sl-2B Z ) 2 . 



(30) 
(31) 

(32) 
(33) 



Then, inserting eq. (|27p into cq. (|26p. one obtains the 
following more explicit expression for the time-optimal 
evolution of the unitary operator, i.e. 



u OPT (t) 



cos ™Ch 3 (t)-sin^B$S%(t) 

i^B S}$(t)[cos<t>(t)al - sin0(t)cr 2 
Sit 



where 4>{t) = ~[6(t) + 0(O)]/2. 

The expression ([Ml) for UoPT(t) still depends upon the 
integration constants £?o, B z , SI and 0(0) and the coupling 
J. These constants, together with the optimal duration 
time T of the evolution and the irrelevant global phase 
X, can be finally fixed by imposing the target condition 

CD). 

We are interested here in the time-optimal realization 
of the symmetric entangler gate Uf 3 (see eq. (4) in [43] ) , 
i.e. 



U f = U : 



13 



(35) 



This is diagonal in the 1, 3 qubit subspace and explicitly 
reads 



U f = e l TUh\ Uh 3 :=Diag(-l, 1,1,1). (36) 



Imposing the target condition eq. (fT0|) with U opt if) 
given by eq. (fM]) and /7/ given by eqs. ([51?]) - ([TO]) , and 
separately equating the terms multiplying the Pauli op- 
erators u 2 ,a 2 and <j\ and the identity operator l 2 , re- 
spectively, we obtain the following set of conditions for 
operators acting on the 1,3 qubit subspace: 



i[x+f]tfl3 



B o cos0(T)^ 3 (T), 
B sm4,(T)S^(T), 
SIT 



(37) 
(38) 



or 



sin ^Ch 3 (T) + cos —B$Sg(Tft<i) 
cos^^(T)-sin^^5L 3 (Tl.40) 



From eq. (|40|) (imposing the reality of its left hand 
side term) we immediately obtain the value of the global 
phase x = (k — l/4)n, where k e Z. Then, the only non 
trivial solution [50] for eqs. ([3"Tl) - ([3"8]) is easily seen to be 



S$(T) = 0, which is possible, using eq. (|28[) , if and only 
if S-)-(T) = s-(T) = sq(T) = and, upon comparing 
with eq. ([30]) - ([51 1 . provided that 



lo±T = im± 1 ujqT = nno, 



(41) 



where n± and no are positive integers. Substituting the 
solution ([4lT) into eqs. ([39]) we also obtain the time- 
optimal integral of the motion SIT = 2irm, where m € Z. 
Finally, inserting the time-optimal formulas for Si and 
w± , wq into eq. (|40p , we obtain the conditions 



k + 2p; 



+ 2r-l; n = n+ + 2q - (L42) 



with fco := fc + m + 1, p,q,r G Z, and the parity of n-o 
and n_ the same and the opposite of that of n + . 

At this point we note that the formulas ([3"2"[) - ([3"3"[) can 
be easily inverted (using ([IT]) ) to write down the param- 
eters Bq,B z and J and the optimal time duration T as 
functions of the integers n± and no. Introducing the fol- 
lowing (odd-integer- valued) functions of the integers n± 
and no in order to simplify the notation, 



/ ± (n + ,n_,n ) := n\ 
/o(n+,n_) := n + 

from eqs. ([32[) - ([3"3"]) we thus obtain 



n 2 _ ± 2n 2 



0) 



2 2 

n , — n_ 



(JT) 2 
(i? T) 2 

iB z T) 2 



f- 



fo 



V2 f 



(43) 
(44) 



(45) 
(46) 

(47) 



We immediately see that one can minimize JT from eq. 
(j45[) by minimizing /_ as a function of the integers k,p, q 
and r (via eqs. ([42]) and ([43]) ). In other words, we 
can determine the optimal time duration T of the quan- 
tum evolution necessary to realize the target gate U[ 3 
as measured in terms of the coupling J. Then, substi- 
tuting the time-optimal values of k,p,q and r into eqs. 
([35 ]) - ([H ]) and subsequently into (|g5 ]) -(|57 ]) . we obtain the 
time-optimal values of the integrals of the motion Bq and 
B z . Finally, we still have to impose the energy constraint 
(j!4p which, in terms of the functions f± and /o, explicitly 
reads 



(wT) 2 



7T 

T 



/+- 



i 



f- 



\/2 /o 
4 



(.48) 



Substituting for the time-optimal values of k,p,q and r 
into eqs. ([4*3"]) - (|4"4"]) and (|35]l, this constrains the values 
of the coupling J as a function of the (given) energy 
parameter uj. 

More explicitly, let us look first for the minima of T 
via eq. (|45p . For a start, we note that, due to the parity 
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properties of and no (eqs. (|42|0 . the function 

/_ (n + , n_, n ) is always an odd integer. Moreover, its 
minimum value /_min = 1 in principle may be achieved 
in two cases, i.e. either when n+ is even (ko even) or 
odd (ko odd). It is immediate to check that the case of 
n+ even (and therefore, via (|42j) . n~ and no odd) is not 
possible [5l|. Furthermore, in the case in which n + is 
odd, defining ko := 2ko + 1 (with ko € Z), the condition 
/— MIN = 1 translates into (p+ko)(2r — iq+l) = 2q 2 — r 2 . 
Substituting /_min = 1 into eq. (|46p we can rewrite 
(B T) 2 = -167T 2 {[q+r+2(p+k )] 2 -l/8}[(q-r) 2 -l/8]. 
Imposing the reality condition, i.e. that (BoT) 2 > 0, 
after a simple algebra we finally obtain that the only 
possible solution for /-min = 1 would be given by (1 — 
2q)[l + 2q + 4(p + ko)] = 1, with r = q or r = -q/(l - 
2q) and q ^ 0,1, which no set of integers {ko,P, q,r} 
can ever satisfy. The next step is to check whether the 
next-to minimum value of /_nm = 3 may be consistent 
with the reality condition for (BoT) 2 > 0. Again, it is 
immediate to verify that the case of n + odd is impossible 
[5^ . However, in the case of n + even, by defining fco := 
2k E (with k E € Z), the condition /_nm = 3 translates 
into (p + k E )(2r - Aq + 1) = 1 - r(r - 1) + 2q(q - 1). 
Calculating (|46l) for /_nm = 3 and n + even, we obtain 
(B T) 2 = -(n 2 /3){[A(p + k E ) + 2q - l] 2 - 3/2}[(2g - 
l) 2 - 3/2]. The reality condition (B T) 2 > now can be 
satisfied either for the values of the integers p + kg = 1, 
with q = 0, —1 and r = 0, —1, or for p + k E = — 1, with 
(7 = 1,2 and r = 1,2. Summarizing, we have found that 
the optimal duration time required to realize a U( 3 gate 
with the Ising Hamiltonian (|11[) is given by 



(49) 



the corresponding time-optimal magnetic fields are 
V5 



\B, 



2^3 
f2_ 



ttJ. 

7tJ 



(50) 
(51) 



while the time-optimal precession frequency of the mag- 
netic field is given by 



mi 



\ I - 7rm J, 



(52) 



where the integer m and the phase 9(0) are still arbitrary. 
Moreover, the coupling J and the energy parameter lu are 
constrained (from eq. (0S])) by the following condition 



J 2V3 



11 + 8 m 



-1/2 



(53) 



IV. MODIFIED HAMILTONIAN AND 
TIME-OPTIMAL CNOT(l,3) 

An analysis similar to that performed in the last Sec- 
tion can be done in the following two situations. On 
the one hand, one can assume to be in the physical sit- 
uation in which the Hamiltonian (|11[) is available to a 
given experimentalist O working with the standard com- 
putational basis, which in our three-qubit model cor- 
responds to the Hilbert space spanned by the states 
{|0), |1)} g> {|0), |1)} <g> {|0), |1)}. Then, we can think of 
another experimentalist O' who performs measurements 
in the basis which corresponds the Hilbert space spanned 
by the states {|0), |1)} ® {|0), |1)} ® {|+), |-)}, i.e. the 
basis where the rotated states {|+), |— )} (with |+) := 
W\0) = (|0) + |1»/V2; I-) := W\l) = (|0> - \l))/y/2 
and W is the Walsh-Hadamard transform) are used for 
qubit 3, while qubits 1 and 2 use the standard computa- 
tional basis. More formally, this situation is equivalent 
to the second experimentalist seeing an effective, rotated 
Hamiltonian given by 

H'(t):= ^[J 12 (t)ala 2 z +J 23 (t)a 2 z a 3 x ]+B(t)-a 2 . (54) 

On the other hand, one may think of the situation in 
which one experimentalist who can perform measure- 
ments in the standard computational basis for the three- 
qubit system, but with the Hamiltonian (|54[) available. 

In both cases, the Hamiltonians (fTTj) and (|54[) are re- 
lated by the transformation 



H'(t) = VH(t)V, 



(55) 



where we have introduced the operator V :— 1®1(8 
W. Then, one may formulate a QB problem using an 
action principle similar to that of the previous Section. 
In other words, one obtains the quantum brachistochronc 
eq. ([9]), and the form of the available Hamiltonian (|54|) 
is guaranteed by the operator F' formally given by eq. 
(|15p with the same set of Lagrange multipliers, but with 
the constraint (|13[) now replaced by 



f 2 := Tr(Fa^)-47rJ = 



(56) 



(while the energy constraint (fT4|) is the same) . A lengthy 
but simple algebra shows that the quantum brachis- 
tochrone equations are again given by (|16p. with the ex- 
ception that now we have to replace everywhere 

Piz -> Pix ; Kiz -> ^zix ; Vi e {x,y, z}. (57) 

Then, we can follow the same procedure as in the pre- 
vious Section and find that eqs. (fTB]) admit the integral 
of the motion B 2 + B 2 = (B' ) 2 and the general solution 
(|17p. In particular, the time-optimal Hamiltonian can be 
seen to become 



H' OPT (t) = VH OPT (t)V, 



(58) 
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where Hopr{t) is given by cq. (fT9)l . The correspond- 
ing time-optimal evolution operator is then given by 
U' OPT {t) = VU PT(t)V, with the Uoprit) of cq. (|34|) . 
and with the operators Bjy, Sjj and C^ 3 always given by 
eqs. ([25]) and p8 ]) -([29 ]l . respectively. 

Now the goal is to time-optimally synthesize the gate 
CNOT(l,3), i.e. 

U f := CNOT(l,3) = e -'t( 1 + CT >x-^-^) j (59) 

which can be diagonalizcd (in the 1, 3 qubit subspace) as 

U f = VU'^V, (60) 

where we have introduced the operator (acting in the 1,3 
qubit subspace) U' D := Diag(l, 1, 1, — 1). Following the 
same methods of the previous Section, we then find that 
the target condition (|10l) is equivalent to impose again 
the operator conditions (f3"T)) - f4"0")) . with the only differ- 
ence that in (|4"0"j) we have to replace exp(ix + 7r/4)/7^ 3 — > 
exp(ix')C^'i) ■ This implies that the optimal value of the 
global phase is given now by x' = k'ft, with k' £ Z, while 
the optimal f2 is the same as in the previous Section. 
The non trivial solution for eqs. (|3"T|) - (J3"5| is again given 
by s+(T) = s~(T) = Sq(T) = 0, and consequently, by 
eq. (|4"Tj) for the same functions uj± and uiq given by eqs. 
(f32"|) - (|33"j) . However, substitution of the optimal values of 
X and f2 and ui± and uiq into eq. (|40p now gives, instead 
of eqs. (0, 

n' + = k' Q + 2p; n'_ = n + + 2r' + 1; n' = n' + + 2g'(61) 

with k' := k' + 771, p' , q' , r' £ Z, and the parity of n' and 
n+ the same and the opposite of that of n'_ . Considera- 
tions similar to those made in the previous Section lead to 
the expressions (|45[) ~ P5)) . and to the minimization of the 
evolution time T at /_nm = 3 (with f± and fo always de- 
fined by (|43l) - (|44)) ). In this case the condition /_nm = 3 
translates into n' + (2r' - Aq' + 1) = 3 - (2r' + l) 2 + 8q' 2 , 
which cannot be satisfied if n' + is even [53[. When n' + 
is odd instead, defining k' Q := 2k' Q (k' £ Z), the con- 
dition /— nm = 3 becomes (p' + k' )(2r' — 4q' + 1) = 
2q'(q' + 1) - r'(r' + 2). Then we can rewrite (B' T) 2 = 
-(167r 2 /3){[2(p' + k' ) + q' + l] 2 - 3/8}(q' 2 - 3/8), and 
imposing the condition (B T) 2 > finally gives q' = or 
q' = -[l+2(p'+k' )] and (2r'+l)[4(p'+k' )+2r'+3] = 3. 
The latter can be satisfied either if p' + k' Q = 0, together 
with r' = 0, -2, or if p'+k' = -1, together with r' = ±1. 

In conclusion, the time-optimal duration necessary to 
realize a CNOT(l,3) gate with the Hamiltonian (|5"4")) is 
given again by (|49p . while it is easy to check that the 
time-optimal magnetic fields |Bo|,|B^|, the precession 
frequency f2 and the ratio J/\oj\ between the coupling 
J and the energy parameter oj are also still expressed by, 
respectively, eqs. 0, (EJ) (EH) and (j53"|) . 



V. DISCUSSION 

Wc have presented the exact and analytical solution for 
the time-optimal realization of two entangling gates, the 
U[ 3 and the CNOT(l, 3), between two indirectly coupled 
qubits, labelled 1 and 3, in a 3-qubit linear spin chain 
subject to, respectively, an Ising type interaction or a 
slightly modified Ising Hamiltonian, where in both cases 
a local magnetic field can be applied on the intermedi- 
ary qubit 2 and the constraint of a finite available energy 
is imposed. No constraints ensuring instantaneous local 
unitary operations are imposed. In particular, we showed 
that the U( 3 gate can be optimally realized via an I sing 
Hamiltonian of the same form as that discussed in [42( 
and that the time required is shorter than that found 
for the particular decomposition of the unitary evolution 
considered in (42[. We then presented the analytical so- 
lution for the time-optimal realization of the CNOT(l, 3) 
via the slightly modified Ising Hamiltonian (j54)) . which 
again is shown to require the same time duration as the 
U{ 3 gate. 

Of course the orbits generated by the Ising Hamilto- 
nian (fTTj) cannot directly reach the target CNOT(l,3), 
which can instead be reached via the modified Ising 
Hamiltonian (|5^|) . According to the standardparadigm 
of time-optimal quantum computing (see, e.g., [2|), where 
one-qubit unitary operations are assumed to have zero 
time cost, one might infer that the CNOT(l,3) can be 
still generated via the Ising Hamiltonian (|11[) provided, 
for example, that a series of additional local operations 
acting on qubits 1, 2 and 3 are applied in the sequence 
given by formula (13) of 0. Somewhat surprisingly, this 
would appear to take the same time ([4T)]) , which is shorter 
than that found in [42j ■ However, the reason of the better 
performance of our quantum evolution stems from the 
fact that the orbit described in [42[ is based on a spe- 
cial decomposition of the unitary operator as a product 
of unitary factors, each of which is indeed time-optimal, 
while the whole product would be time-optimal only if 
these single factors commuted which each other, which is 
not true in general. 

On the other hand, in the more physical and general 
QB framework, where the time cost of local unitarics 
is not negligible in principle, one can instead think of 
the new Hamiltonian (fM]) as generated by a change of 
basis of the Hilbert space for the three-qubit system ac- 
cording to {|0>, |1>}® {|0>, {|0>, |1)} -> {|0>, 

{|0), |1)} ® {|+), |— )}, as explained in the previous Sec- 
tion. This would relate the Ising Hamiltonians (TTTI) and 
(|S"4")) via cq. (f551) . In other words, one can now see that 
the time-optimal realization of the gate CNOT(l,3) un- 
der the modified Hamiltonian (p>4"| and where both the 
control and the target qubits belong to the standard com- 
putational basis {|0), |1)} is, in fact, equivalent to the 
time-optimal realization of the gate CNOT ± (l,3) under 
the Ising Hamiltonian (fTTj) where the control qubits be- 
long to the basis {|0), |1)} while the target qubits belong 
to the basis {|+), |-)}. 
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We should note that our problem is slightly more re- 
stricted than that of [42] , in the sense that we also impose 
the finite energy condition eq. ((6]). In fact, the normal- 
ization condition Q constrains not only the amplitude 
\B\ of the local magnetic field in the Hamiltonian (fTTj) 
(or ([54]) ). via cqs. ([50]) and (f5Tj) . but also the allowed 
values for the Ising coupling constant J, via eq. (|53|) . 
As a consequence, only a (discrete) set of possible values 
for the coupling constant J are actually allowed, and the 
optimal value of J is related to the energy available ui 
through eq. (|53[) . In particular, one can see that, in the 
case when the energy uj available increases (for a fixed 
integer m) the coupling J grows and, therefore, via eq. 
(|4"9"|) , the optimal time required to generate the gates de- 
creases, while the strengths of the optimal local magnetic 
field, via eqs. ([50[) and (f5Tj). and also the precession fre- 
quency, via eq. (|52[) . grow (which is what is physically 
expected). Alternatively one can use eq. ([53)) to express 
the optimal values of |i?o| and \B Z \ (eqs. ([50)) and ([51]) ) 
as a function of the total energy available uj, and then 
to evaluate the ratio between the amplitude of the local 
magnetic field controlling the time-optimal evolution and 
the value of the Ising coupling constant J available in the 
experiment, obtaining 



\B\ 



J 2V3\ 



5 + 8 m 



(62) 



The ratio (|62[) is a monotone increasing function of the 
integer m, and it tells us that the amplitude of the 
driving magnetic field cannot be smaller than |B|min = 
2vr(l - v /378) 1 / 2 /(\/3J) (for m = 1). Furthermore, al- 
though the fact that the allowed values of J belong to a 
discrete set (determined by the values of the integer m 
via eq. ([53])) might appear as a limitation of our model, 
we point out that this feature is just the reflection of our 
simplifying technical choice (made just for ease of pre- 
sentation) of considering only equality constraints (J5J. 
We conjecture that this feature should be mitigated by 
extending the QB methods to the more general and phys- 
ical ansatz where, e.g., the total energy available is only 
bounded from above (i.e., substituting the equality con- 
straint ([5] with an inequality constraint, and using the 
proper variational techniques, see, e.g., [54J). 



Furthermore, we would like to comment on the relation 
of our results with those of [3l|. In particular, in [3l[ the 
one-qubit components of the time optimal Hamiltonian 
were found to be constant, which is not the case in the 
present work. However, in j3lj there is no constraint for 
the one-qubit or the two-qubit components of the Hamil- 
tonian, except for the normalization, which results in the 
absence of such components in F' (eq. ([5])) and thus 
yields the constancy of the one-qubit components. Here 
we have the constraints (|12p ([TB"]) on the one-qubit and 
the two-qubit components of the Hamiltonian and there- 
fore the result changes. 

Our work was in part motivated by [42j , and therefore 
we specialized our analysis to the case of the Hamilto- 
nians ([TT]) and (fM[) . where the local operations are in 
principle allowed only on the intermediary qubit. How- 
ever, there is no particular reason why one should limit 
to such a case, and one may think instead of more gen- 
eral situations in which also local interactions for (one 
of) the boundary qubits are available. Similarly, one 
could extend the analysis to the case of different interac- 
tion couplings between the indirectly coupled qubits (i.e. 
Ji2 J23, see [HI) or for different topological couplings 
among the qubits (see, e.g., [l?} and [HI). Similarly, a 
simple and straightforward generalization of our results 
would be to include the study of how to time-optimally 
generate quantum gates other than the CNOT (e.g., the 
SWAP(1, 3), the ,/SWAP(l,3) etc.). A further natural 
extension of our research will be also to consider our QB 
methods in the study of linear (Ising) chains of n qubits 
(see (3)), e.g. in the search for time-efficient ways of cre- 
ating and propagating spin order and coherences along 
the chains, or the transfer speed of single spin excitations 
via external magnetic fields along Heisenberg spin chains 
(see, e.g., [55[ and references therein). 
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